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ABSTRACT 

A  method  is  presented  to  compute  the  span  wise  vorticity  in  polar  coordinates  from  2D  vertical  cross  sections 
of  high-resolution  line-of-sight  Doppler  wind  lidar  observations.  The  method  uses  the  continuity  equation  to 
derive  the  velocity  component  perpendicular  to  the  observed  line-of-sight  velocity,  which  then  yields  the 
spanwise  vorticity  component.  The  results  of  the  method  are  tested  using  a  ground-based  Doppler  lidar,  which 
was  deployed  during  the  Terrain-Induced  Rotor  Experiment  (T-REX).  The  resulting  fields  can  be  used  to 
identify  and  quantify  the  strength  and  size  of  vortices,  such  as  those  associated  with  atmospheric  rotors. 
Furthermore,  they  may  serve  to  investigate  the  dynamics  and  evolution  of  vortices  and  to  evaluate  numerical 
simulations.  A  demonstration  of  the  method  and  comparison  with  high-resolution  numerical  simulations 
reveals  that  the  derived  vorticity  can  explain  66%  of  the  mean-square  vorticity  fluctuations,  has  a  reasonably 
skillful  magnitude,  exhibits  no  significant  bias,  and  is  in  qualitative  agreement  with  model-derived  vorticity. 


1.  Introduction 

The  Terrain-Induced  Rotor  Experiment  (T-REX) 
field  campaign  was  carried  out  in  spring  2006  in  Owens 
Valley,  California  (Grubisic  et  al.  2008).  The  main  fo¬ 
cus  of  the  experiment  was  the  investigation  of  the 
coupled  mountain  wave  rotor  system  in  the  lee  of  the 
Sierra  Nevada.  One  of  the  primary  objectives  is  to 
study  the  internal  dynamics  of  rotors  and  their  inter¬ 
action  with  boundary  layer  processes  near  the  earth’s 
surface. 

Atmospheric  rotors  are  characterized  by  strong  down- 
slope  surface  winds  in  the  lee  of  a  mountain  range  that 
decelerate  rapidly  to  a  weaker,  transient,  and  turbulent 
flow,  often  directed  back  toward  the  mountain  (e.g., 
Holmboe  and  Klieforth  1957;  Kuettner  1959;  Doyle 
and  Durran  2002;  Doyle  et  al.  2009).  Based  on  3D  high- 
resolution  numerical  simulations,  Doyle  and  Durran 
(2007)  proposed  that  the  largest  hazard  for  aviation  is 
actually  not  the  whole  rotor  but  smaller  vortices,  which 
they  called  subrotors.  These  vortices  can  be  produced 
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along  a  shear  line  that  is  generated  by  boundary  layer 
separation  resulting  from  strong  downslope  winds  and 
lifted  aloft  by  the  lee  wave  circulation  (Doyle  et  al. 
2009).  Some  of  the  major  challenges  for  T-REX  are  to 
observe  these  subrotor  vortices,  to  estimate  the  strength 
of  the  horizontal  vorticity,  and  to  investigate  their 
characteristics. 

In  this  study,  high-resolution  scanning  Doppler  lidar 
observations  taken  during  T-REX  are  analyzed,  and 
a  method  for  deriving  tangential  velocity  and  vor¬ 
ticity  fields  from  lidar  line-of-sighf  velocities  Vr  is  pre¬ 
sented  (see  Fig.  1  for  the  definition  of  the  coordinate 
system).  During  T-REX  rotor  events,  most  lidar  obser¬ 
vations  were  made  in  a  vertical  cross-sectional  mode 
along  a  baseline  260°  from  north  to  80°  from  north, 
which  is  roughly  parallel  to  the  main  flow  and  perpen¬ 
dicular  to  the  valley  axis  and  the  Sierra  Nevada  Ridge. 
The  observations  were  made  by  the  Institut  fiir  Physik 
der  Atmosphare  of  the  Deutsches  Zentrum  fiir  Luft- 
und  Raumfahrt  (DLR)  using  a  scanning  2-/rm  coherent 
Doppler  lidar  system  from  Lockheed  Martin  Coherent 
Technologies  (WindTracer;  details  available  online  at 
http;//www.lockheedmartin.com/ssc/coherent/products/ 
windtracer/Speciflcations.html).  The  lidar  was  situated 
on  the  alluvial  slope  of  the  Owens  Valley  at  36.79°N, 
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Fig.  1.  Schematic  diagram  of  the  Cartesian  and  cylindrical  co¬ 
ordinate  systems.  The  x-z  plane  and  the  r-<f>  plane  are  chosen 
perpendicular  to  the  valley  axis  (i.e.,  the  x  axis  is  pointing  toward  an 
azimuth  of  80°  from  north,  and  the  y  axis  is  pointing  toward  350° 
from  north).  Lidar  scans  were  performed  mainly  in  the  x-z  plane 
and  the  r-<^  plane,  which  are  roughly  parallel  to  the  main  flow. 


118.21°W,  and  1240  m  MSL,  which  is  approximately 
2  km  west  of  the  town  of  Independence.  The  valley  is 
confined  by  the  Sierra  Nevada  to  the  west,  which  in¬ 
cludes  several  peaks  above  4000  m  MSL  and  the  highest 
peak  in  the  lower  United  States,  and  the  Inyo  Mountains 
to  the  east,  with  peaks  exceeding  3000  m  MSL.  The 
width  of  the  valley  is  approximately  20-25  km.  The  time 
required  for  a  full  lidar  cross  section  was  between  30 
and  150  s  during  the  T-REX  field  campaign,  depending 
on  scanning  speed.  The  cross  sections  used  in  this  note 
were  taken  over  a  time  period  of  90  s.  The  resolution  of 
this  cross  section  is  105  m  along  the  lidar  beam  and 
2°  elevation. 

Doppler  lidars  only  measure  the  wind  component 
in  the  direction  of  the  lidar  beam,  whereas  either  u  and 
V  or  Vy  and  are  needed  to  derive  spanwise  vortic- 
ity  fields.  Under  relatively  stationary  flow  conditions, 
this  can  be  solved  by  observations  with  two  (or  more) 
ground-based  lidar  systems  or  by  making  observa¬ 
tions  from  several  directions  with  a  lidar  mounted  on 
a  moving  platform  (e.g.,  an  aircraft;  Reitebuch  et  al. 
2003).  Atmospheric  rotors  during  T-REX,  however, 
were  highly  transient  and  associated  with  strong  tur¬ 
bulence. 

In  addition  to  the  DLR  lidar,  Arizona  State  Uni¬ 
versity  (ASU)  operated  a  similar  Doppler  lidar  during 
the  T-REX  field  campaign.  Both  lidars  performed  co- 
planar  lidar  observations  with  the  goal  of  deriving 
the  horizontal  and  vertical  velocity  through  a  combi¬ 
nation  of  both  datasets  through  dual  Doppler  analysis 
(Drechsel  et  al.  2009;  Hill  et  al.  2010).  Results  from 
dual  Doppler  can  also  be  used  to  derive  spanwise 
vorticity,  but  the  combination  of  two  different  datasets 
inevitably  leads  to  errors  resulting  from  the  time  dif¬ 
ference  of  the  observations  (up  to  about  30-150  s  with 
the  scan  strategy  during  T-REX).  Furthermore,  the 
volume  of  the  observations  is  never  completely  the 


same  because  of  the  observation  geometry.  Thus,  it  is 
concluded  that,  in  events  with  strong  turbulence  and 
transient  flow,  the  time  sequencing  and  differences  of 
the  observational  volume  likely  lead  to  significantly 
larger  errors  than  the  method  presented  in  this  note. 
For  example,  relatively  small  vortices  propagate  with 
the  mean  flow,  which  can  lead  to  a  displacement  of  up 
to  600  m  (nearly  six  range  gates)  in  30  s  for  a  wind 
speed  of  20  m  s^\ 

This  note  presents  an  alternative  method  that  exploits 
the  continuity  equation  to  derive  and  the  spanwise 
vorticity  from  vertical  slice  line-of-sight  observations  of 
a  single  lidar  (section  2).  This  approach  has  the  advan¬ 
tage  that  the  time  difference  between  two  adjacent  ob¬ 
servations  is  in  the  range  of  only  0.2-1  s.  The  errors  of 
lidar  V,.  observations  are  small  compared  to  the  velocity 
fluctuations  associated  with  turbulent  flows  as,  for  ex¬ 
ample,  in  rotors.  Errors  of  coherent  2-/rm  Doppler  lidars 
are  often  estimated  to  be  on  the  order  of  0.1  m  s^^  un¬ 
der  favorable  conditions  and  may  increase  by  a  factor 
of  approximately  2-3  at  a  lower  signal-to-noise  ratio 
(Grund  et  al.  2001;  Chai  et  al.  2004). 

The  main  errors  for  the  derived  and  vorticity  fields 
are  expected  to  occur  because  of  flow  divergence  per¬ 
pendicular  to  the  plane  of  the  lidar  cross  section.  For  this 
reason,  the  errors  of  the  calculation  resulting  from  di¬ 
vergent  lateral  flow  are  investigated  with  tests  using  2D 
and  3D  high-resolution  numerical  simulations  with  the 
atmospheric  portion  of  the  Coupled  Ocean-Atmosphere 
Mesoscale  Prediction  System  (COAMPS)  in  section  3. 
Furthermore,  the  application  of  the  calculation  to  T-REX 
observations  is  also  shown  in  section  3.  A  discussion  of 
the  results  is  presented  in  section  4. 


2.  Tangential  velocity  and  vorticity  calcnlations 
derived  from  Vy  observations 


The  tangential  velocity  (i.e.,  velocity  component 
perpendicular  to  Vy)  can  be  calculated  for  the  grid  points 
in  each  quadrant  through  the  use  of  the  2D  continuity 
equation  in  polar  coordinates  and  the  application  of 
a  boundary  condition  such  that  the  vertical  velocity  is 
zero  at  the  surface.  It  could  also  be  calculated  for  the 
whole  hemicycle  starting  at  (^  =  0°,  but  errors  of  the 
calculation  accumulate  with  every  elevation  increment. 
It  follows  that  it  is  more  accurate  to  begin  the  calculation 
from  both  sides  (i.e.,  at  cj)  =  0°  and  cj)  =  180°,  re¬ 
spectively). 

Assuming  that  there  is  no  divergence  perpendicular  to 
the  plane  of  the  lidar  observations  (r-cf)  plane),  the 
continuity  equation  in  polar  coordinates  (r,  cj))  is 


1  5 
r  dr 


(rVy)  + 


r  flch 


=  0, 


(1) 
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where  r  is  the  radial  distance  from  the  lidar;  0  is  the 
elevation  angle;  and  Vr  and  are  the  radial  and  tan¬ 
gential  velocities,  respectively.  At  the  point  (/,  k),  dV^ldr 
can  be  approximated  as 

dr  2Ar  ’  ^ 

where  Ar  is  the  radial  distance  between  two  range  gates, 
i  is  the  index  of  the  range  gate,  and  k  is  the  index  of 
the  elevation  step.  Similarly,  dV^ldcf)  at  the  point  (i,  k) 
can  be  approximated  as 

^  ~  ^KQ', <:-!)] 

ad)  Ad)  ■  ^ 


Using  approximations  (2)  and  (3),  Eq.  (1)  can  be  mod¬ 
ified  to 


r(i,k) 


rAd) 

2A7 


r(i+\,k) 


-  V 


r(i-l,k)J 


(4) 


Equation  (4)  can  be  solved  stepwise  by  starting  at 
(f)  =  0°,  where  _  i)  =  0,  and  incrementing  to  0  =  90°. 
Likewise,  can  be  calculated  from  cf)  =  180°  to  </>  =  90°. 
However,  in  the  case  of  divergence  in  the  y  direction, 
this  leads  to  errors  of  dV^ldr  resulting  in  errors  of  that 
propagate  from  one  elevation  increment  to  the  next. 
This  aliasing  effect  causes  “rings”  of  high/low  and 
large  gradients  of  dV^/dcf).  To  mitigate  the  aliasing, 
a  smoothing  parameter  S  is  introduced  to  average  at 
the  grid  point  (i,  k  —  1),  resulting  in 


rAO 

-  Ad)U^^.^^  —  2^[^r(/+l,A:)  “ 

(5) 

Tests  with  model  and  lidar  fields  show  that  the  best 
results  are  achieved  with  S  in  the  range  of  0.2-0.33  (see 
section  3).  No  significant  improvements  are  achieved 
with  a  smoothing  algorithm  that  includes  the  second 
level  of  range  gates  (not  shown).  This  smoothing  algo¬ 
rithm  is  expected  to  reduce  errors  that  arise  as  a  result  of 
the  numerical  approximations  in  Eqs.  (2)  and  (3)  and  as 
a  result  of  the  small-scale  variability  of  the  flow.  How¬ 
ever,  the  smoothing  procedure  does  not  remove  errors 
arising  from  large-scale  lateral  flow  divergence. 

Without  the  along-beam  smoothing,  the  method 
produces  spurious  structures  of  in  some  locations,  as 
mentioned  previously.  This  is  understandable,  because 
errors  resulting  from  lateral  divergence  in  one  range 


gate  will  propagate  to  higher  elevations  at  constant  Ar 
because  of  the  iterative  algorithm.  However,  divergence 
resulting  from  small-scale  flow  variability  is  usually 
compensated  by  convergence  nearby.  Thus,  the  spurious 
structures  can  be  effectively  mitigated  through  the 
smoothing  procedure,  which  is  confirmed  by  the  visual 
inspection  of  the  results  that  show  realistic  structures 
and  no  obvious  rings  after  smoothing.  The  improve¬ 
ments  through  such  an  approach  are  further  supported 
by  the  tests  with  output  fields  of  the  numerical  model. 

As  a  consequence  of  larger-scale  lateral  flow  diver¬ 
gence,  there  is  a  difference  of  at  0  =  90°  for  the  cal¬ 
culations  beginning  with  0  =  0°  and  0  =  180°.  Therefore, 
a  correction  is  performed  by  assuming  a  constant  error  of 
[K(/+i,A:)  “  ^r('-iA)]’  which  is  equal  to  the  difference  of 
the  two  calculations  at  0  =  90°  divided  by  the  number 
of  elevation  steps  (further  discussed  in  section  3). 

After  calculating  the  vorticity  field  can  be  corn- 
puled  by  using  the  definition  of  the  spanwise  vorticity 
component  tj  in  polar  coordinates, 


1  5 
rSr 


(6) 


For  the  vorticity  calculation,  dV^ldr  and  dVrld(f>  at  the 
point  (/,  k)  were  approximated  using  centered  finite 
differencing,  similar  to  that  in  Eq.  (2). 


3.  Evaluation  of  the  velocity  and  vorticity 
calculations  using  model  fields  and 
T-REX  observations 

An  example  of  T-REX  lidar  observations  is  shown 
in  Fig.  2.  A  more  detailed  discussion  of  the  rotor  event 
on  16  April  2006  can  be  found  in  Doyle  et  al.  (2009).  The 
observations  were  taken  during  one  of  the  strongest 
rotor  events  of  T-REX  and  show  strong  downslope 
winds  about  8-10  km  upstream  of  the  lidar.  A  region 
with  weaker  flow  partly  directed  backward  toward  the 
mountain  range  can  be  seen  between  the  lidar  and  5  km 
upstream  of  the  lidar.  Strong  westerlies  are  apparent 
above  this  region  with  weaker  flow.  A  distinct  vortex 
with  a  diameter  of  about  1  km  that  formed  at  the  shear 
layer  at  about  3  km  MSL  is  visible  2  km  upstream  of  the 
lidar. 

A  vertical  slice  from  the  3D  high-resolution  numerical 
simulations  using  an  eddy-resolving  version  of  the 
CO  AMPS  model  was  used  to  optimize  the  smoothing 
parameter  S  to  decide  if  the  along-beam  smoothing  and 
the  correction  for  differences  at  0  =  90°  should  be  ap¬ 
plied  and  to  estimate  the  uncertainty  of  the  and 
vorticity  calculation.  The  grid  spacing  of  the  simulations 
is  60  m  in  the  horizontal  and  50  m  in  the  vertical  (see 
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f)  LOS  velocity  model 


a)  LOS  velocity  lidar 


b)  Tangential  vel.  lidar  without  correction 


g)  Tangential  velocity  model 
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c)  Tangential  velocity  lidar  with  corr. 
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d)  Vorticity  lidar  without  correction 
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e)  Vorticity  lidar  with  correction 
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h)  Tan.  velocity  calculated  using  continuity 
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i)  Vort.  in  polar  coord,  with  model  vtan 
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j)  Vort.  in  polar  coord,  with  vtan  using  cont. 
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Fig.  2.  (a)  Here  Vy  is  measured  by  the  DLR  Doppler  lidar  during  T-REX  at  2134  UTC  16  Apr  2006.  The  scan  is  a  vertical  cross  section 
along  a  baseline  from  (left)  260°  to  (right)  80°,  which  is  roughly  perpendicular  to  the  valley  axis.  Observations  were  taken  during  an  event 
with  strong  downslope  winds  to  the  lee  of  the  Sierra  Nevada.  Cold  colors  indicate  flow  toward  the  lidar,  and  warm  colors  indicate  flow 
away  from  the  lidar.  The  topography  is  shown  with  a  solid  black  line,  (b)  Here  is  calculated  from  (a)  using  the  continuity  equation,  S  = 
0.33,  and  no  correction  for  differences  at  </>  =  90°;  (c)  as  in  (b),  but  calculated  with  correction  for  differences  at  <^  =  90°;  (d)  vorticity  using 
(a),(b);  (e)  vorticity  using  (a),(c);  (f)  Vy  from  model  simulation  along  the  same  line  as  (a);  (g)  from  same  model  simulation  as  (f);  (h) 
calculated  from  (f)  using  the  continuity  equation,  S  =  0.33,  and  a  correction  for  differences  at  </>  =  90°;  (i)  vorticity  calculated  using  (f),(g); 
and  (k)  vorticity  calculated  using  (f),(h).  (All  velocities  are  in  m  s^*,  and  all  vorticity  values  in  s^^.) 
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Table  1.  Standard  deviation  and  mean  of  vorticity  and  fields. 
The  calculations  were  made  with  the  correction  for  differences  at 
<j>  =  90°.  (The  units  of  velocity  are  m  s^^,  and  the  units  of  vorticity 
are  s~^.) 


Std  dev 

Mean 

Vorticity  from  model  interpolated 

0.0298 

0.0058 

Vorticity  using  model  tangential  velocity 

0.0278 

0.0059 

Vorticity  using  continuity,  S  =  0.2 

0.0270 

0.0070 

Vorticity  using  continuity  S  =  0.33 

0.0256 

0.0070 

V,!,  from  model 

15.74 

11.62 

V,f,  using  continuity,  S  =  0.2 

15.22 

13.30 

using  continuity,  S  =  0.33 

15.01 

13.30 

Doyle  et  al.  2009).  The  baseline  of  the  cross  section  is 
the  same  as  for  the  most  common  scan  pattern  during 
T-REX  (roughly  parallel  to  the  main  flow).  There  are 
obviously  differences  between  the  simulations  (Fig.  2f) 
and  the  lidar  observations  (Fig.  2a),  but  the  simulations 
appear  to  be  a  suitable  test  case  because  the  strength  and 
size  of  the  vortices  have  a  similar  magnitude  and  struc¬ 
ture.  More  details  of  the  simulations,  similarities,  and 
differences  to  T-REX  observations  are  discussed  in 
Doyle  et  al.  (2009). 

The  model  output  fields  were  averaged  to  a  polar  grid 
that  is  similar  to  the  lidar  observations  with  a  radial 
resolution  of  105  m  and  2.0054°  for  0.  In  the  vicinity  of 
the  lidar,  where  the  resolution  of  the  polar  grid  is  higher 
than  that  the  model  grid,  the  angular  tolerance  for  the 
averaging  was  increased  by  a  factor  of  4  for  the  first  four 
range  gates  and  by  a  factor  of  2  for  range  gates  6-10. 
After  averaging  the  model  winds  to  a  polar  grid,  Vr  and 
were  calculated  (Figs.  2f,g).  Additionally,  was 
calculated  from  model  Vr  using  the  continuity  equation, 
as  discussed  in  section  2  (F,^-CONT;  Fig.  2h). 

The  spanwise  vorticity  was  calculated  from  a  vertical 
slice  of  the  model  test  fields  using  three  different 
methods:  1)  Cartesian  coordinates  and  then  averaged  to 
polar  grid  (not  shown  in  Fig.  2);  2)  polar  coordinates 
using  Eq.  (6)  with  approximations  described  by  Eqs.  (7) 
and  (8)  (Fig.  2i);  and  3)  polar  coordinates  using  V,j,  cal¬ 
culated  with  the  continuity  equation  (Fig.  2j;  VORT- 
CONT).  The  calculations  for  Figs.  2h  were  made  using 
S  =  0.33  with  a  correction  for  differences  of  V^2lI^  =  90°. 

The  standard  deviation  and  mean  of  these  fields 
are  shown  in  Table  1.  All  values  of  Tables  1,  2,  and  3 
were  only  calculated  for  the  domain  where  corrected 


F^-CONT  is  available  (Fig.  2j)  to  make  the  values 
comparable.  Table  2  shows  the  mean-square  differences 
of  calculated  using  the  continuity  equation  to  the 
averaged  model  for  a  range  of  S  from  zero  (no 
smoothing)  to  0.33.  The  setting  of  5  =  0.33  implies  that 
the  three  range  gates  used  for  the  averaging  are  nearly 
equally  weighted.  In  addition,  these  fields  are  normal¬ 
ized  with  the  mean-square  values  of  Table  1.  Table  3 
shows  the  vorticity  differences  in  a  similar  manner  as 
the  differences  displayed  in  Table  2. 

The  and  vorticity  fields  are  calculated  with  sim¬ 
plified  numerical  approximations  of  the  derivatives, 
which  inevitably  leads  to  numerical  errors  proportional 
to  the  grid  increment.  For  this  reason,  the  vorticity  cal¬ 
culated  in  polar  coordinates  has  smaller  fluctuations 
than  the  vorticity  calculated  in  Cartesian  coordinates 
(Table  1).  The  calculation  of  using  continuity  and  the 
smoothing  parameter  S  leads  to  a  further  decrease  of 
fluctuations  that  can  be  resolved  by  the  method.  FIow- 
ever,  small  differences  of  the  vorticity  fields  calculated 
in  Cartesian  coordinates  and  the  ones  in  polar  co¬ 
ordinates  (not  shown)  indicate  that  these  numerical  er¬ 
rors  are  smaller  than  the  ones  resulting  from  divergence 
perpendicular  to  the  lidar  scan. 

Figure  2h  shows  that  the  calculation  of  V,/,  using 
continuity  leads  to  reasonably  accurate  results  in  the 
region  with  fairly  uniform  flow  (e.g.,  upstream  3  km  and 
above  4  km  MSL),  but  that  lateral  divergence  can  lead 
to  errors  of  for  example,  in  the  turbulent  region 
above  the  lidar,  where  the  model  fields  show  reversed 
flow  in  contrast  to  calculated  F^.  Despite  these  errors  of 
F^-CONT,  the  VORT-CONT  in  the  same  area  is  in 
reasonable  agreement  with  the  vorticity  fields. 

Overall,  the  error  of  F^(,-CONT  is  6%  of  the  total 
fluctuations  for  S  >  0.1  with  the  correction  for  differ¬ 
ences  at  0  =  90°  and  7%  for  S  &  0.2  without  correction 
(Table  2).  Both  smoothing  and  the  correction  show 
a  clear  improvement  of  F^-CONT.  According  to  the 
results,  S  should  be  in  the  range  of  0.2-0.33  with  lowest 
errors  at  5  =  0.33.  The  total  fluctuations  of  F^-CONT 
(Table  1)  are  realistic  given  the  fact  that  the  test  field  is 
highly  turbulent. 

The  smoothing  is  assumed  to  reduce  errors  that  pri¬ 
marily  arise  because  of  divergence  perpendicular  to  the 
plane  of  the  lidar  scans  associated  with  turbulent  fluc¬ 
tuations.  As  a  consequence,  the  optimal  setting  of  S  is 


Table  2.  Differences  of  calculated  from  model  Vr  using  the  continuity  equation  with  and  without  a  correction  for  differences  3X<j>  =  90° 

and  model  V^.  The  values  are  normalized  with  values  shown  in  Table  1. 


Smoothing  parameter  S 

0.00 

0.10 

0.20 

0.30 

0.33 

Relative  mean-square  difference  (F^  with  correction) 

0.08 

0.06 

0.06 

0.06 

0.06 

Relative  mean-square  difference  (F^  without  correction) 

0.10 

0.08 

0.07 

0.07 

0.07 
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Table  3.  Differences  of  vorticity  calculated  using  the  continuity  equation  with  and  without  a  correction  for  differences  at  4>  =  90°  and 

vorticity  calculated  from  model  velocities  averaged  to  polar  grid. 


Smoothing  parameter  S 

Mean-square  difference  (vorticity  with  correction)  in  10^“*  m^ 
Relative  mean-square  difference  (with  correction) 

Mean-square  difference  (vorticity  without  correction)  in  10^"^  m^  s”^ 
Relative  mean-square  difference  (without  correction) 


0.00 

0.10 

0.20 

0.30 

0.33 

5.82 

3.60 

3.01 

2.78 

2.73 

0.72 

0.45 

0.37 

0.34 

0.34 

7.32 

4.00 

3.22 

2.90 

2.85 

0.91 

0.49 

0.40 

0.36 

0.35 

expected  to  be  dependent  on  the  ratio  of  the  dominant 
turbulence  length  scale  and  the  resolution  of  the  ob¬ 
servations.  The  optimal  S  parameter  is  defined  as  the 
value  of  S  that  leads  to  the  smallest  square  difference  of 
U^-CONT  and  VORT-CONT  to  model  and  vor¬ 
ticity,  respectively.  This  dependency  on  S  is  demon¬ 
strated  for  the  3D  simulation  fields  shown  in  Fig.  3, 
where  the  optimal  S  parameter  is  shown  for  different 
along-beam  resolutions. 

Further  tests  using  a  2D  simulation  (Fig.  3)  show  that 
the  smoothing  parameter  also  improves  the  results  in  the 
absence  of  divergence  perpendicular  to  the  cross  sec¬ 
tion.  Otherwise,  the  optimal  S  parameter  would  be  zero 
or  close  to  zero  for  2D  model  fields.  This  is  interpreted 
as  a  reduction  of  numerical  errors  by  the  smoothing, 
which  also  explains  why  the  optimal  S  parameter  shows 
a  clearly  lower  dependency  on  the  along-beam  resolu¬ 
tion  for  the  2D  fields.  A  direct  comparison  of  the  errors 
using  2D  and  3D  output  fields  is  not  possible  because  the 
flow  field  itself  is  significantly  different  in  2D  and  3D 
simulations. 

The  error  resulting  from  the  vorticity  calculation  is 
less  than  37%  for  S  &  0.2  with  the  correction  and  less 
than  40%  without  correction.  The  smallest  error  of 
VORT-CONT  is  reached  using  S  =  0.33.  With  this 
setting,  the  error  of  the  calculation  is  34% ,  which  implies 
that  the  vorticity  derived  using  the  continuity  equation 
can  explain  66%  of  the  vorticity  fluctuations.  The  posi¬ 
tive  effect  of  the  correction  actually  does  not  appear  to 
be  that  large  for  the  vorticity  field  compared  to  the  im¬ 
provement  of  V^-CONT.  This  is  understandable,  be¬ 
cause  only  the  gradients  and  not  the  absolute  values  of 
are  important  for  the  vorticity.  A  qualitative  com¬ 
parison  of  the  results  shown  in  Fig.  2  indicates  reason¬ 
able  agreement  of  the  vorticity  fields:  the  magnitude  of 
the  fluctuations  of  VORT-CONT  is  realistic  (Table  1), 
and  the  maxima  are  of  similar  magnitude  (—0.2  s^^). 

The  application  of  the  method  to  T-REX  observations 
is  shown  in  Figs.  2a-e.  The  qualitative  interpretation  of 
using  different  smoothing  parameters  for  vorticity  cal¬ 
culations  with  T-REX  observations  confirmed  that  0.2  < 
S  <  0.33  is  a  reasonable  choice  for  the  parameter  (not 
shown).  Because  the  calculation  of  needs  a  continu¬ 
ous  field,  five  isolated  missing  values  of  were  in¬ 


terpolated  from  neighboring  values  for  the  calculation 
of  and  vorticity. 

Despite  differences  of  that  arise  as  a  result  of  the 
correction,  the  results  for  vorticity  are  very  similar  with 
and  without  the  correction  for  differences  at  0  =  90° 
(Fig.  2).  Because  observations  at  all  elevation  angles  are 
needed  to  perform  a  correction,  the  coverage  of  cor¬ 
rected  VORT-CONT  is  lower  than  that  of  the  un¬ 
corrected  field.  In  the  field  displayed  in  Fig.  2,  this  is 
clearly  a  limitation,  because  the  strongest  vortex  ap¬ 
pearing  at  a  height  of  2  km  AGE  1.7  km  upstream  of  the 
lidar  is  not  completely  covered  in  the  corrected  field. 
The  vorticity  calculation  has  been  tested  with  several  of 
the  lidar  scans  during  T-REX  and  shows  reasonable 
results  in  all  cases.  These  tests  confirm  that  the  cal¬ 
culation  leads  to  reasonable  values  for  vorticity.  The 
V,^-CONT  field  is  improved  further  with  the  inclusion 
of  the  correction  factor. 

4.  Discussion 

We  have  demonstrated  that  realistic  vorticity  fields 
can  be  derived  from  vertical  slice  lidar  Vr  scans  with 
a  relatively  simple  and  computationally  inexpensive 
method  that  calculates  solving  the  2D  continuity 


Along-beam  resolution  (m) 

Fig.  3.  The  S  parameter  leading  to  the  lowest  errors  of  the  and 
vorticity  calculation  as  a  function  of  along-beam  resolution  in  tests 
with  2D  (solid  lines)  and  3D  (dotted  lines)  model  output  fields. 
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equation  numerically.  An  along-beam  smoothing  of 
is  essential  to  retain  realistic  results.  A  further  correc¬ 
tion  for  differences  of  at  cj)  =  90°  with  a  constant  error 
term  improves  the  fields  of  but  does  not  appear  es¬ 
sential  for  the  calculation  of  vorticity  because  only  er¬ 
rors  of  the  along-beam  gradient  of  lead  to  vorticity 
errors  and  not  errors  of  the  absolute  V^. 

Tests  with  an  eddy-resolving  model  output  from  a  3D 
simulation  reveal  that  the  vorticity  derived  using  the 
continuity  equation  can  explain  66%  of  the  vorticity 
fluctuations.  In  addition,  the  fluctuations  of  vorticity 
have  a  realistic  magnitude,  and  the  calculation  does  not 
show  a  significant  bias.  The  derived  vorticity  and  model 
vorticity  are  qualitatively  in  good  agreement.  The 
field  calculated  from  using  the  continuity  equation 
has  an  accuracy  of  94%,  but  errors  are  larger  in  regions 
of  strong  turbulence. 

The  method  is  based  on  the  assumption  of  2D  conti¬ 
nuity,  whereas  the  flow  also  has  fluctuations  in  the  third 
dimension  perpendicular  to  the  vertical  slice.  Thus,  the 
main  errors  in  the  derived  tangential  velocity  and  vorticity 
result  from  small-scale  flow  variability  and  large-scale 
divergence  in  the  third  dimension.  The  first  contribut¬ 
ing  factor  should  be  mitigated  through  the  along-beam 
smoothing.  This  is  shown  by  tests  with  the  velocity 
fields  from  the  model  simulations  where  the  third  ve¬ 
locity  component  is  known.  The  errors  resulting  from 
large-scale  flow  divergence  are  attempted  to  be  reduced 
through  a  correction  based  on  the  differences  of 
at  0  =  90°  after  calculating  the  tangential  velocity  in¬ 
crementally  from  the  surface  on  both  sides  of  the  lidar 
up  to  0  =  90°.  Without  a  large-scale  flow  divergence,  the 
calculations  from  both  sides  should  lead  to  the  same 
tangential  velocity  at  cf)  =  90°.  Thus,  the  field  was 
linearly  relaxed  to  minimize  the  difference  with  this 
correction.  Some  of  the  remaining  differences  in  the 
field  indicate  that  the  correction  does  not  fully  remove 
these  errors.  However,  as  illustrated  by  the  model  sim¬ 
ulations,  this  large-scale  flow  divergence  primarily  af¬ 
fects  the  tangential  velocity,  and  the  derived  vorticity  is 
not  substantially  impacted.  Instrumental  noise  associ¬ 
ated  with  coherent  Doppler  lidars  is  estimated  to  be  in  a 
range  of  0.1-0.3  m  s^^,  according  to  Grund  et  al.  (2001) 
and  Chai  et  al.  (2004);  thus,  it  is  comparably  smaller  than 
the  velocity  gradients  of  several  meters  per  second  be¬ 
tween  neighboring  range  gates.  Errors  resulting  from 
experimental  noise  can  also  be  expected  to  be  reduced 
by  the  along-beam  smoothing,  because  these  errors  are 
spatially  uncorrelated. 

In  summary,  the  proposed  method  is  useful  to  quan¬ 
tify  the  strength  and  size  of  vortices,  such  as  those  as¬ 
sociated  with  atmospheric  rotors,  which  pose  a  severe 
hazard  to  aviation.  One  limitation  of  the  method  is  that 


only  the  spanwise  vorticity  component  is  derived.  How¬ 
ever,  because  there  is  no  observing  system  that  can  ob¬ 
serve  the  full  3D  flow  field  within  a  time  range  of  minutes, 
it  is  seen  as  a  useful  approach  for  topographic  flow  prob¬ 
lems  and  other  flows  when  the  main  wind  component  is 
along  the  plane  of  vertical  slice  lidar  scans. 

The  along-beam  smoothing  is  performed  through  a 
weighted  moving  average  of  the  neighboring  range  gates, 
whereby  different  weighting  parameters  for  the  average 
are  investigated.  The  optimal  smoothing  parameter  is 
expected  to  be  dependent  on  eddy  size,  resolution  of  the 
observations,  and  accuracy  of  the  observing  system.  This 
smoothing  may  need  to  be  adapted  when  the  method  is 
applied  for  different  observing  systems  or  if  the  vortices 
of  interest  are  significantly  different,  but  the  general 
approach  can  be  assumed  to  be  valid  for  various  appli¬ 
cations  where  the  strength  of  vorticity  needs  to  be 
quantified  by  vertical  cross  sections  of  line-of-sight  ob¬ 
servations  from  Doppler  lidar  or  also  Doppler  radar 
instruments.  In  recent  years,  several  airports  around  the 
world  (e.g.,  Hong  Kong  International  Airport;  Shun  and 
Chan  2008)  installed  operational  wind  lidar  systems  for 
aircraft  safety,  and  the  presented  method  is  also  seen  as 
a  possibility  to  provide  quantitative  information  on 
vortices  and  potential  aviation  hazards  in  real  time. 
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